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Abstract 

The exact renormalization group equation for pure quantum gravity is used 
to derive the non-perturbative /3-functions for the dimensionless Newton 
constant and cosmological constant on the theory space spanned by the 
Einstein-Hilbert truncation. The resulting coupled differential equations are 
evaluated for a sharp cutoff function. The features of these flow equations 
are compared to those found when using a smooth cutoff. The system of 
equations with sharp cutoff is then solved numerically, deriving the complete 
renormalization group flow of the Einstein-Hilbert truncation in d = 4. The 
resulting renormalization group trajectories are classified and their physical 
relevance is discussed. The non-trivial fixed point which, if present in the 
exact theory, might render Quantum Einstein Gravity nonperturbatively 
renormalizable is investigated for various spacetime dimensionalities. 



I. INTRODUCTION 



Exact renormalization group (RG) equations [Q] provide a powerful tool for the non- 
perturbative investigation of both fundamental (renormalizable) and effective quantum 
field theories. In particular the RG equation of the effective average action |^ has been 
applied to a variety of matter field theories as well as to Euclidean quantum gravity . 

The effective average action is a Wilsonian (coarse grained) free energy functional 
defined on a smooth spacetime manifold. Its construction is based upon a modified version 
of the standard path integral for the generating functional which has a built-in infrared 
(IR) cutoff at the mass scale k. All quantum fluctuations with momenta > k'^ are 
integrated out as usual, while the contributions from modes with p"^ < k'^ are suppressed 
by the cutoff. This cutoff is implemented by giving a momentum dependent (squared) 
mass Rkip^) to the modes with covariant momentum p [Q]. 

The fc-dependence of is governed by an exact functional RG equation. In any real- 
istic theory it is impossible to solve this equation exactly, but by appropriately truncating 
the space of action functionals ("theory space") one can obtain nonperturbative approx- 
imate solutions which do not rely upon small expansion parameters. The truncation is 
carried out by making an ansatz for which contains finitely or infinitely many free, k- 
dependent parameters gi{k). Upon inserting this ansatz into the functional RG equation 
and projecting the RG flow onto the truncation subspace one obtains a coupled system 
of ordinary differential equations for the gj. 

For the case of Euclidean quantum gravity, the effective average action and its RG 
equation have been constructed in where also the flow equations for the "Einstein- 
Hilbert truncation" were derived. The Einstein-Hilbert truncation retains only the in- 
variants / d'^Xy/gR and / d'^'x^/g, the associated coupling constants being Gk, the running 
Newton constant, and A^, the running cosmological constant, respectively. The original 
construction of ref. employs a cutoff of "TYPE A" which is formulated in terms of 
the complete metric fluctuation /i^,^. In refs. [^,0 a new cutoff of "TYPE B" was in- 
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troduced which is formulated in terms of the component fields of hf^,^ appearing in its 
transverse-traceless decomposition |0. 

One of the interesting predictions of the Einstein-Hilbert truncation is that the high 
energy behavior of 4-dimensional quantum gravity is governed by a non-trivial ultravio- 
let (UV)-attractive fixed point for the dimensionless Newton constant and cosmological 
constant ^J^J^,^]. (In d dimensions, they are defined as Qk = k'^~'^Gk and = Afc/fc^, 
respectively.) If this fixed point is present also in the exact theory, Quantum Einstein 
Gravity is likely to be renormalizable at the nonperturbative level 0. Despite its per- 
turbative nonrenormalizability it could then be a predictive, fundamental theory valid 
at arbitrarily small distances. An UV fixed point g* for entails that Gk = g*/k'^~'^ 
vanishes for k —>■ oo (provided c? > 2) so that the theory becomes asymptotically free. 

Also the "phenomenological" implications of the running gravitational constants for 
black hole physics and for cosmology |]12|,ll3l have been investigated. In particular. 



it has been argued [0 that the UV fixed point reflects itself in the cosmology of the 
Planck era and that it might lead to a solution of the flatness and the horizon problems 
of classical Friedmann-Robertson- Walker cosmology. 

Previous investigations and applications of the Einstein-Hilbert truncated flow equa- 
tions have mainly been based on approximate solutions of the RG-equations in the far 
IR {k 0) and far UV {k oo). Besides these solutions only little is known about the 
complete RG flow given by the Einstein-Hilbert truncation. 

In this paper we shall therefore study the numerical solutions to the flow equation for 
the dimensionless cosmological constant and the dimensionless Newton constant gk, 
analyzing the RG flow in the complete X-g— plane. For a generic cutoff, the pertinent f3- 
functions /3;^ and I3g contain rather complicated "threshold functions" which are functions 
of Afc and gk, and functionals of Rk{p'^)- In order to make a numerical solution of the flow 
equation feasible we introduce a sharp cutoff ior which the integrals deflning the 

threshold functions can be evaluated analytically. In contrary to many standard matter 
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field theories where the singular nature of the sharp cutoff often leads to ill defined or 
divergent /3-functions, we shall see that applying the sharp cutoff to the /3-functions arising 
form the Einstein-Hilbert truncation leads not only to well defined functions f3y^,(3g, but 
in addition yields the same characteristic behavior of the RG trajectories. Moreover, for 
universal quantities, even the quantitative results are very similar to those found when 
using a smooth cutoff function. 

The /3-functions employed in the present investigation are those derived in the original 
paper 0, i.e. they arise from a cutoff of "TYPE A" and a constant gauge fixing parameter 
a = 1. 

The remaining sections of this paper are organized as follows. In Section II we give a 
brief summary of the original derivation 0] of the "TYPE A" flow equation for and Gk-, 
introducing all the definitions needed for the later analysis. In Section III we investigate 
the properties of the threshold functions as an important ingredient of the flow equation. 
In this course we introduce a sharp cutoff function which allows the analytical evaluation 
of the integrals appearing in these functions. The Sections IV and V contain the results 
of our numerical analysis. In Section IV we investigate the renormalization group flow 
below the Planck scale and compare the results obtained with the sharp and a smooth 
exponential cutoff function. In Section V we use the sharp cutoff to derive the full RG 
flow on the (7-A-phase space and classify all possible solutions to the flow equation in 
d = 4. Here we also investigate the reliability of the Einstein-Hilbert truncation in ci = 4 
and its limitations for > 4 by investigating the behavior of the non-trivial fixed point 
in various other dimensions. 

II. THE FLOW EQUATION OF THE EINSTEIN-HILBERT TRUNCATION 

In order to derive the non-perturbative flow equation for the dynamical Newton con- 
stant Gk and the cosmological constant on the theory space spanned by the Einstein- 
Hilbert truncation we use the effective average action approach to pure quantum gravity. 
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For the details of the following discussion we refer to the original paper P]. 

The main ingredient of this method is the exact evolution equation for the effective 
average action rfc[(7^j,] for gravity. The derivation of this evolution equation parallels the 
approach already successfully tested for Yang-Mills theories . In principle it is also 

possible to include the additional renormalization effects of and Gk coming from matter 
fields | jr7| , p!B|] , but these are not included in the present derivation. 

In the construction of Tk[g\ one starts out with the usual path integral of d-dimensional 
Euclidean gravity. This is gauge fixed by using the background field method [ p!9| , p0[ 
and employing a background gauge fixing condition. Due to the introduction of the 
background metric g the effective average action Tk[g]g\ now depends on both the full 
metric g and the background metric g. The conventional effective action T[g] is regained 
as the /c — > limit of Tk[g] = '^k[g] 9 = g], where the two metrics have been identified. In 
this manner Tk[g] becomes invariant under general coordinate transformations. 

The crucial new component in the construction of rk[g,g] is the fc-dependent IR- 
cutoff term AkS in the action under the path integral. This term discriminates between 
the high {p^ > k"^) and low-momentum modes {p^ < k"^). It suppresses the contribution of 
the low-momentum modes to the path integral by adding a momentum dependent mass 
term 

AkS[h, C, C; -g] = ^k' J d'x ^h,,Rr''\gY''"'K'^ + d'^x ^ C,Rt[g]C^ (2.1) 

Here is a constant, and the first and second term on the RHS provide the cutoff for 
the fluctuations of the metric h^i, = g^y — g^u and the ghost fields C^, C^, respectively. 
In this work we choose the following form of the cutoff operators Rf^^ and Rf^: 

Rrm = zrk'R^'\-D/e), Rt[g] = eR^'\-D/e) (2.2) 

Here (Zf''")^^''" = [(/ - P^)'"^^" - (rf - 2)/2P^'^''"] Z^k is a matrix acting on /i^^. In this 
expression (P^)'^'^^" = d^^ g^^ g^" projects h^^ onto its trace part (p. In the terminology of 
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this form of the cutoff function defines the "cutoff of TYPE A". The so-called shape 
function is arbitrary except that it has to satisfy the conditions 



2; ^ oo 







(2.3) 



Neglecting the evolution of the ghost sector which corresponds to a first truncation of 
the general structure of Ffc, one finds that rk[g,g] satisfies the following flow equation: 



dt^k [g, 



It. 
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-Tr 



>^-'T^^^ + nr [g]r'd,Rr [9] 



M[g,g]+Rtm 'd^RTlg] 



(2.4) 



(2) 

Here T^. [g, g] denotes the Hessian of Tk [g, g] with respect to g^i, at fixed background field 
(jfj^u and t = ln{k/k) is the "renormalization group time" with respect to the reference 
scale k. Furthermore, Ai represents the Faddeev-Popov ghost operator. 

In order to obtain the nonperturbative flow equation for the running Newton constant 
Gk and the cosmological constant we now approximate Tk[g, g] by the Einstein- Hilbert 
truncation, considering the subspace spanned by the operators / d'^Xy/g and / d'^Xy/g R 
only: 



Tk[g,g] =il6nGk)-' Jd'^xy/^ { -R + 2Afc| + classical gauge fixing 



(2.5) 



Substituting this ansatz into the evolution equation (|2.4| ) and projecting the resulting flow 
onto the subspace given by the Einstein-Hilbert truncation then leads to the flow equation 
for Afc and Gk- This equation is conveniently written down using the dimensionless Newton 
constant g^ and cosmological constant A^: 



For these couplings the flow equation reads 



Afe = \k ^ 



(2.6) 



dt\ = ^xi>^k, gk), dtgk = /3„(Afe, gk) 



(2.7) 
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where the /3-functions are given by 



(3,iX,g) = -{2-r)^)X + -{i7:y-'^/'g 



2 d{d + 1) $i/2(-2A) - 8 d $i/2(0) - d{d + 1) $^/2(-2A) 



I3gi\,g) = {d-2 + riN) g 

Here 77 at is the anomalous dimension of the operator / d'^x^R, 

gB,{\) 

and the functions -Bi(A) and i?2(A) have the following definition: 



(2.8) 



(2.9) 



d{d + 1) $Vi(-2^) - - 1) *^/2(-2A) 
-4d$i/2_i(0)-24<|.^/2(0) 

S2(A) ^ (47r)^-'^/2 + 1) *i/2_i(-2A) - 6d(d - l)*^/2(-2A) 



(2.10) 



Furthermore, we introduced the i?*^"-* -dependent threshold functions and (p 



1,2, 



as 



T{n) Jo 



1 



[z + i?(o)(z)+H'' 



n— 1 



i?W(^) 



(2.11) 



r(n)yo [z + R(o){z)+wY 

In order to investigate the renormalization group flow below the Planck scale it is more 
convenient to introduce a second set of coupling constants which are made dimensionless 
by using an arbitrary but fixed scale M rather than k: 



G(k) = Gk M"-^ , X(k) = XkM 



-2 



(2.12) 



For convenience we also introduce k = kM~ as a dimensionless scale parameter. 
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Substituting the dimensionless coupling constants ( p.l2| ) into the /3-functions 
finds the following set of differential equations 

dk 



8|) one 



(2.13) 



k -txCk) =VnX + 1 (47r)i-'^/2 ^ 2 d{d + 1) ^^^/^{-2X/P' 
dk 2 



8rf<l>^/2(0)-rf(rf+l)r7^<l>^/2(-2AA 



d/2\ 



with 



Vn 



k'^-^GBiiX/k^ 



(2.14) 



(2.15) 



1 - k'^-^GB2{\/k'^) 

Here dt = kd/dk = kd/dk has been used. The equations ( p.l3|) , ( p. 141) and ( p. 15]) will be 
our starting point for investigating the renormalization group flow below the scale M in 
Section IV. 



III. THRESHOLD FUNCTIONS 



Before we can turn to the numerical solution of the flow equations ( |2.8| ) and ( p.l3| ), 
( p.l4| ) it is important to understand the characteristics of the threshold functions $^ and 
$^ which depend on the form of the cutoff function R^^^ chosen. 



A. Smooth Cutoff Functions 



For the purpose of studying the general properties of the threshold functions, we 
flrst assume that satisfles (g^ and is a smooth function which does not vanish too 
quickly for small values of its argument: z + R^^\z) > 1, V^; > 0. Under this condition 
the deflnition of the threshold functions (|2.11|) immediately shows that they vanish for 



w 



oo. 
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Furthermore one notes by examining the denominator of the integrands in (|2.11|) that 
the functions and are well defined only for arguments w > —1. For w < —1 

the denominator vanishes for some z in the region of integration, yielding a non-integrable 
singularity. Therefore and are finite and well defined for w G (—1, oo) only. 

One can derive a recursion formula which connects the threshold functions with dif- 
ferent p-values. Interchanging the derivative with respect to w and the 2;-integration 
yields: 

= -pK^\w), ^^li^) = -pK^\^) (3.1) 



B. The Exponential Cutoff 

For practical computations it is necessary to have an explicit form of the cutoff function 
R^^\ In the calculations done in |^,|ri],|T2| the one-parameter family of exponential cutoffs 

i?(°)(z;s)^"P = for s>0 (3.2) 

exp(s z) — 1 

has been used. Here s is a "shape parameter" . For these cutoffs the threshold functions 
( [^.11| ) are seen to be positive definite: 

<l>P(u;)^"P>0, $^(w)^"P>0, p= 1,2,... , n = 1,2,... and w e (-1,00) (3.3) 

The integral defining the threshold functions with exponential cutoff can be carried out 
analytically for vanishing argument. Using the integral representations for polylogarithms 



and the Riemann ^-function, Li,^(s) = /q°° dz and ({u) = ItT '^^ 7^ lH 
one easily verifies 

*n(0)'^^' = 7^{C(^+l)-LWl(l-^)} 



^liOf'-'' = —^7^ rLi„_i(l - .) (3.4) 



1 

But for non-vanishing arguments w an analytic solution to these integrals is unknown 
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C. The Sharp Cutoff 



In order to be able to evaluate the threshold functions for any argument w we introduce 
a different cutoff function, the sharp cutoff. On the level of the dimensionful function 
Rkip^) = k"^ R^^^p^ / k"^) the sharp cutoff is defined as 

R,{p^r^RQ(l-p^/e) (3.5) 

where the limit i? ^ oo is to be taken after the integration over p, i.e. after substituting 
( p.5|) into the threshold functions. In the path integral this choice of cutoff leads to a 
complete suppression of modes with momentum p^ < k"^ while all modes with p"^ > k"^ are 
completely integrated out without the additional mass term. 

The evaluation of the integrals in the threshold functions then proceeds as follows. In 
the first step one substitutes z = p^ /k"^ and w = vjkP' into the definition of 

V[n) Jo [p^ + Rk[p^) + vY 

Here the prime denotes the derivative of R^^^ with respect to its argument. The integrand 
of equation ( |3.6| ) can then be written as a total derivative with respect to /c, for any p > 1: 



2T{n){p-l)^ ' Jo D k {p^ + Rk{p^) + v)P-^ 

(3.7) 



^ {ky--^rdp^k-^ 



2r(n)(p-l)' ' Jo dk {p^ + Rk{p^)+v)P-^ 

Here the derivative D / Dk acts by definition only on the fc-dependence of Rk-, but not on 
V. In order to rewrite it in terms of d/dk we introduced the constant v which is strictly 
independent of k. Only at the end of the calculation it is identified with v = wk'^. (For 
p = 1 the formula (|3.7| ) breaks down since its right hand side is no longer well defined.) 
Assuming p > 1, we interchange the 2;-integration and the derivative with respect to k. 



This is allowed since the integral in ( |3.7| ) is absolutely convergent. If one now substitutes 
the sharp cutoff Rk{p'^Y^ = -R6(l —p^jk^) one finds: 
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^>/kT = -Wwd^T-^ k— dp 



2T{n){p-l) dk Jo {p^ + Re{l-pyk^)+v)p-i 
Taking the limit R ^ oo restricts the momentum integration to p"^ < k"^: 



(3.8) 



(3.9) 



2r(n)(p-l) aA; Jfca (p2 + {))p-i 

The resulting integral is trivially evaluated by acting with the /c-derivative on the lower 
integration limit. This yields our final result for the threshold functions with a sharp 
cutoff: 

mwY'' = -^^—- ^— — for p>l (3.10) 

' r(n) p-l {l + w)P~^ ^ ^ ' 



One easily verifies that these threshold functions also satisfy the recursion relation ( |3.1| ) 
derived for a smooth cutoff. This relation can be used to define (^^iwy^ as the solution to 
the differential equation -^^^iwY^ = y+w "which arises from substituting $^(w)'^'^ 

into (|3.1|). ^l^i^wy^ is then determined up to a constant of integration: 

= - 1^(1 + ^) + ^ni^r (3.11) 

r(n) 

The constants $^(0)'^'^ = ipn are undefined a priori; they parameterize the residual cutoff 
scheme dependence which is still present after having opted for a sharp cutoff. 

For numerical calculations we need to fix the V9„'s. Usually we shall assume them 
equal to the corresponding constants arising from the exponential cutoff with the shape 
parameter s = 1: 

^ KiOr = Kiof"' ^'^'^ (3.12) 

This choice leads to a very good agreement between the threshold functions evaluated 
by using the exponential and the sharp cutoff, at least for values of w in the "confidence 
interval" [—0.7,1]. This is shown in FIG. |I|, where the numerical values of those 
functions which appear in the /3-functions for d = 4 are compared for the two types of 
cutoffs. 
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Exp s=l 




Exp s=l 




FIG. 1. Comparison between the threshold functions in the /3- functions for d = A for exponential 
and the sharp cutoff. The horizontal lines indicate the "confidence interval" in which the two cutoff types 
yield similar results. 

In order to investigate the residual cutoff scheme dependence of the results obtained 
from the sharp cutoff it is useful to define a one-parameter family of v?n's, analogous to the 
one-parameter family of exponential cutoff functions ( p.2D . This is done by generalizing 
the relation ( p.12 ) to arbitrary values of the shape parameter s: 

^„(s)^$i(0f^P(^) (3.13) 

The resulting functions v'i(s) and ^2{s)i which are the only constants of integration 
appearing in d = 4, are shown in FIG. |^. Here one finds that for s > 15 the numerical 
values of (pi{s) and (p2{s) are not subject to relevant changes any more. Therefore we will 
limit further investigations of the residual cutoff scheme dependence to the region s < 15. 
For the $^-functions the integral appearing in ( p. 11 ) can also be evaluated explicitly 
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s=30 



25 30 



FIG. 2. Dependence of the constants of integration (pi and ip2 on the shape parameter s. 

by substituting the sharp cutoff. In a calculation similar to the one shown above one 
obtains 



for p > 1 



for p = 1 



(3.14) 



Having at hand the explicit form of the threshold functions, equations ( p.lO| ), ( |3.11| ) 
and (|3.14|) , we are now in a position to write down the /3-functions ( p.8|) with all integra- 
tions carried out. The result reads 



^g{\,gr = id-2 + r^]^)g 

/3,(A, gr = -(2 - Vn) a + \ i^ny-'^' g ■ 



2d{d+l) did+1) 

ln(l - 2A) + 2d{d - 3)ipd/2 - , Vn 



r(rf/2) 



r(rf/2 + i) 



with 



l-gB^iXy 



(3.15) 



(3.16) 



Here i?i(A)'^'^ and B2{XY'^ are given by 
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BliXr = I {^7ry-'/'[ - ^^^^^ -2X) + d{d- 3)^,/2_i 

_6d{d-l) 1 _ 24 1 
T{d/2) 1-2X~ r(rf/2)i 

^^(A)^^ = -i (47r)i-^/2 ^(^ +1) (3.17) 



Thus the flow equation ( p.TD has now boiled down to a coupled system of ordinary differ- 
ential equations whose RHS is known explicitly and which is easily solved numerically. 

IV. RG TRAJECTORIES BELOW THE PLANCK SCALE 



In this section we discuss numerical solutions to the "M-scaled" differential equations 
( ^.13| ), ( p^.l4D for G and A in = 4 dimensions. We specify initial conditions at some 



fixed scale k = k and identify the mass scale M, which was used in ( |2.12|) for defining the 



dimensionless variables G and A, with this initial point M = k. It will prove convenient 
to rewrite the flow equations in terms of the scale parameter 

y = k^ = k^M^ (4.1) 

Hence the initial point k = M corresponds to y = 1. 

We shall see that some RG trajectories can be continued down to A; = 0, while others 
terminate in a singularity at a nonzero value of k. For every trajectory which can be 
extended to = we define a Planck mass in terms of the final, i.e. infrared, value of 
Newton's constant: 

mpi = Go = G{k = 0)-^/' (4.2) 
As a consequence of this definition, 

(5(0) = GqM^ = M^/ml^ (4.3) 
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After having solved the RG equation and having found the final value G(0) one could in 
principle use equation ( [4 .31) in order to express M in terms of the more physical Planck 
mass. 

In order to disentangle the various effects which contribute to the running of G and 
A it is helpful to start the discussion by analyzing two approximate forms of the system 
of equations ( |2.13| ), ( |2.14| ). At a first level of approximation, we neglect the running of G 



by setting ?7jv = and focus on the scale dependence of A alone. The motivation is that 
according to canonical dimensional analysis the running of A is much more "relevant" than 
that of Newton's constant. At a second level of approximation, we further simplify the 
remaining equation, (|2.14|) with rj^ set to zero, by neglecting the backreaction which the 



changing A has on the flow via the threshold functions; to this end we set $^(— 2A/A;^) ~ 

A. The Naive Renormalization Group Flow 

Let us start with the "naive" RG flow which is defined by the following two approxi- 
mations: 

G{k) ^ (5(0) = const, i.e. rj^ 0, and 

^l{~2\{y)/y) ^ (4.4) 
The remaining differential equation = ^ G{fi) is easily solved: 

\{y) = m + $^(0) G(0) (y' - f) (4.5) 

4:71 ^ ^ 

This corresponds to the result of a naive one loop calculation with an IR cutoff which 
also yields a running of A^ proportional to y"^ (x k^. 

The solutions ( [4.5| ) are plotted in FIG. ^ where the dimensionless Newton constant 
was chosen to be G(0) = 1. By equation ([4.3|) this choice implies that M = mpi, i.e. 
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FIG. 3. Solution ( [4.5[) to the naive flow equation for different initial values X{y) and G(0) — 1. 



that the initial conditions are imposed at A; = mpi, and that the cosmological constant is 
measured in Planckian units. 

As FIG. shows, the solutions ( [4.5|) allow for any value of the "renormalized" cosmo- 
logical constant A(0) in the limit k ^ 0, depending on the initial value set at the scale 
y. We could "fine-tune" this initial value so that the renormalized cosmological constant 
becomes A(0) = 0. But this fine-tuning depends on the numerical value of $i(0) which in 
turn depends on the non-physical cutoff function R^^^ chosen. Therefore it cannot have a 
universal meaning. 



B. The Impact of the Threshold Functions 



Now we drop the second assumption made in ( [4. 41 ) and allow for a non-trivial argument 
of the threshold function $2, while keeping rj^ = 0. The resulting flow equation for a 
generic smooth cutoff reads 



^ = ^G{0) [5$^(-2A/i/)-4$^(0) 



For the sharp cutoff it becomes: 

d\ y 



G{0) -5 ln(l - 2A/i/) + (^2 



(4.6) 



(4.7) 



dy 2n 

We shall solve eqs. ( [4.6[ ) and ( [4.7D numerically and compare the resulting trajectories. 
We specify initial values A(y) at the initial point y = y = 1. In order to visualize the 
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FIG. 4. Numerical solutions to the flow equations (4.6) and (4.7) for the exponential and the sharp 
cutoff for G(0) = 1 and shape parameter s = 2. The straight line indicates the boundary X ~ y/2 
introduced by the threshold functions. Trajectories leading to negative values A(0) as well as trajectories 
terminating in the singularity are found. 



generic solution to these flow equations we have to abandon our choice s = 1 for the shape 
parameter because for this special value the RG trajectories have certain properties which 
are not typical. Let us start with s = 2 instead which illustrates the general situation. 
(For the sharp cutoff this change of s leads to a different numerical value of the v'n's, 
whose s-dependence has been defined in equation ( p.l3| ).) 

The resulting trajectories for s = 2 are shown in FIG. ^ The special case arising from 
s = 1 is displayed in FIG. Again we chose (5(0) = 1, implying that M = mpy for those 
trajectories for which G{k = 0) is defined. 

Comparing the "Exp" to the "sc" diagrams in FIGs. ^ and ^ clearly indicates that 
the sharp and the smooth exponential cutoff yield essentially the same RG trajectories. 

The most striking feature of FIG. ^ is that some trajectories cannot be continued 
below a certain finite value y = ytcrm- Trajectories which, in FIG. have led to positive 
IR values A(0) > now terminate because they hit the singular line X = y/2. It is due 
to the singularities of the and ^-functions at w = —2\/y = —1. Looking at the 
flow equations and ( |4.7|) we see that their right hand sides diverge, f3^ +oo, if 
A approaches y/2. Some solutions to eq. ([4.7| ) in the vicinity of A = y/2 are shown in 
FIG. ^ The trajectories approach the singularity with a slope 13^ smaller than 1/2. Only 
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FIG. 5. Behavior of trajectories close to the boundary A = y/2. 

directly on the boundary the slope would jump to oo discontinuously. Equations ( [4.6| ) 
and (|4.7| ) show that the singularity of f3^ also extends to 

A > I (4.8) 

Hence the RG flow is defined only for A < y/2. This leads to the following consequences 
for all cutoffs: 

(i) It is not possible to choose initial values \{y) larger than y/2. 

(ii) There are no trajectories in the region A > y/2. Therefore we obtain only negative 
values for the cosmological constant Aq in the limit k 0. Trajectories that in 
FIG. ^ have led to positive values A(0) > now either lead to negative values of 
the cosmological constant, A(0) < 0, or run into the boundary line X = y/2 and 
terminate at finite values yterm > 0. 

Comparing those trajectories that in FIG. ^ und FIG. ^ yield negative IR-values 
A(0) < we see that including the nontrivial argument of the threshold functions leads 
to a (rather weak) focusing of the cosmological constant towards zero, in the sense that 
the trajectories in FIG. ^ curve upward for y 0. 

Looking at FIG. ^ we see that the choice s = 1 leads to a non-generic behavior since 
the termination of the trajectories at A = y/2 does not occur. In this case all admissible 
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FIG. 6. Numerical solutions to the flow equations (4^) and (4/7) for the exponential and the sharp 
cutoff for G(0) = f and shape parameter s — 1. The bold, straight line indicates the boundary A — y/2 
introduced by the threshold functions. In contrast to the generic case s = 2 no trajectories terminate on 
the boundary A = y/2. 

initial conditions imposed at y = 1 give rise to trajectories which can be continued down 
to y = 0. They all yield a negative or vanishing Aq- 



C. The Complete System 



Let us now also drop the approximation Gk = const and consider the full flow equations 
and ( |2.14D , treating both G{k) and A(A;) dynamical. We compute the solutions for 
the initial values 



G{y = 1) = 0.25, \{y = 1) arbitrary 



(4.9) 



The resulting trajectories for the sharp cutoff are shown in FIG. |^. Here two classes of 
solutions immediately become apparent. 

Trajectories which end at a negative value of the cosmological constant A(0) < have 
already appeared in FIGs. ^ and Later on they will be referred to as of "Type la" . 
Comparing the trajectories with A(0) < of FIG. ^ to those in FIG. ^ and ^ one finds 
that including a dynamical Newton constant in the flow equation counteracts the effect of 
focusing |Ao| towards smaller values, rather leading to even more negative values of A(0). 
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FIG. 7. Solution to the full flow equation with G{y = 1) = 0.25 and various initial values X{y = 1) 
for the sharp cutoff with s = 1. The bold line indicates the boundary at X = y/2. 

The second class of solutions in FIG. ^ which will be classified as Type Ilia, is formed 
by the trajectories ending on the boundary line A = y/2. For A y/2 the anomalous 
dimension r]N rapidly diverges. This leads to a vast increase of the Newton constant G{y) 
preceding the termination of the trajectory. Due to the divergence of 1]^, the RHS of 
both eq. ( p.l3| ) and ( p.l4| ) approaches — oo as A y/2. Therefore the tangents of the 
functions G{y) and \{y) turn vertical causing the trajectory to terminate at a finite value 

2/term > 0. 

Furthermore, FIG. ^ clearly shows the "anti-screening" character of pure quantum 
gravity, i.e. the monotone decrease of the Newton constant Gk with increasing values of 
y or A; 0. 

The most important change arising from the inclusion of a running G{y), which can 
not directly be deduced from the FIGs. ^, || and ^, is the modification of the backward 
evolution when we increase y and try to send it to infinity. For the fiow equations ([4.6|) 



and ( |4.7| ) with G{y) = G{0) kept constant we find that the backward evolution becomes 
undefined for sufficiently large values oiy > y. All trajectories terminate at the boundary 
line A = y/2 at a finite value y < oo and cannot be continued to ?/"="cxd. 

As we shall see in the next section this behavior changes drastically when the running 
of G{y) is included. This is due to the appearance of a non-trivial fixed point which 
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FIG. 8. Solution to the full flow equation with G{y = 1) = 0.25 and various initial values X{y = 1) 
for the exponential with s = 1. The bold line indicates the singularity at A ~ y/2. 



governs the RG flow of tlie coupled G-A— system for large y ^ 1. As a consequence, all 
trajectories shown in FIG. ^ where the running of G was included have a well defined 
backward evolution and can be continued up to ?/"="cx3. 

Next we solve eqs. (|2.13| ) and (|2.14|) using the initial conditions (|4.9| ) and the expo- 
nential cutoff. The results are displayed in FIG. Comparing the trajectories in FIGs. 
1^ and H one sees that the trajectories obtained using the sharp and the exponential cutoff 
are very similar except when they get close to the singular line A = y/2. 

A new phenomenon occurs if we choose initial conditions close to the boundary A = 
y/2. Here the trajectory obtained using the sharp cutoff shows no particularities while the 
one found using the exponential cutoff crosses all the other Type Ilia trajectories which 
start at lower \{y) and terminates at an "unnaturally low" value y, see the first diagram of 
FIG. ^ This behavior is due to the divergence of rj^"^ on a certain 2-dimensional surface 
in A-G-y— space which is discussed in detail in Appendix A. It leads to a termination of 
the trajectories slightly before the A = ?//2-line is reached. This difference between the 
sharp and the smooth cutoff occurs in a region very close to the singularity where the 
Einstein-Hilbert truncation is unreliable in any case. It is clear that it has no physical 
significance. For the sharp cutoff the analogous singularity r/^ — > oo is located in the 
(probably unphysical) region with negative G{y). 



20 



V. THE COMPLETE RENORMALIZATION GROUP FLOW 



After having studied the properties of the RG flow below the scale M in the last section 
we now investigate the complete RG flow, i.e. the RG flow on the entire A-t?— plane. In 
this course we first look at the qualitative features of the flow equation, determining its 
fixed points and their stability properties. Afterwards we construct the full phase portrait 
of the flow by numerically solving the "fc-scaled" flow equation with the sharp cutoff. 

A. The Fixed Points 

1. General Remarks 

The existence and consequences of a non-trivial fixed point in the Einstein-Hilbert 
truncation of pure quantum gravity have already been discussed in refs. Ref. |^ 

uses the exponential shape function ( p.2|) for a cutoff of TYPE A, establishing the existence 
of this fixed point for a wide range of shape parameters s. In Q the properties of the fixed 
point are investigated using the TYPE B cutoff for various smooth shape functions R'^^\ 
In this subsection we extend these surveys to the sharp cutoff within the original 

TYPE A cutoff scheme and compare our results to those obtained by using 

In order to investigate the fixed points we now turn to the flow equations written 
in terms of the /c-scaled coupling constants \,g of (|2.8|) and (|3.15|) . The existence of a 
fixed point A*, (7* requires that the /3-functions of the Einstein-Hilbert truncation vanish 
simultaneously: 

/3,(A = \\g = g*) = 0, /3,(A = X*,g = g*) = (5.1) 

Equation ( |5.1| ) has two solutions: A trivial one for which X* = 0, g* = and a non-trivial 
solution with g* ^ 0. They give rise to a "Gaussian" and a "non-Gaussian" fixed point, 
respectively. 
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An important property of a fixed point is the anomalous dimension r]i\f{X*,g*) at this 
fixed point. From ( pl8| ) and ( pl9D one finds T]^{\*,g*) = and r]N{X*,g*) = 2 - c/ for the 
"Gaussian" and the "non-Gaussian" fixed point, respectively. 

In order to investigate their stability properties we linearize the RG flow at the fixed 
points, 



B ^ [5, 



dX dg 
a\ dg 



(5.2) 



In our case the generalized couplings gi are given by gi = A, g2 = g, and the derivatives 
defining Bij are taken at gj = g*. At an arbitrary point {\,g), and for any R^^\ the 
partial derivatives of the /3-functions are easily found by making use of the recursion 
formula 

dr]N 



dX 



dg 



-(2 - Vn) + (a - |(47r)i-'^/^c/(d+ l)$i/,(-2A)^ 
+1 {4nY-'/' {M {d + 1)^%,{-2X) -2d{d+ 1) m $^/2(-2A 



2 



A_- (47r)i-'^/2d(d+l)$i/2(-2A) 



+ o (4 



,l-d/2 



9(? ■ 2 

{2d [d + l)$i/2(-2A) - 8d $^/2(0) - + 1) r^^v $^/2(-2A) 

:5i'(A)+wi?2'(A)) 



_ g' 
dX l-gB2{X) 



dg 



d-2+\2+ 



9B2{X) 
1-^7 52(A), 



(5.3) 



Here the derivatives of r^jv are given by 



+ 



52(A) 



(i?/(A) + r/^i?2'(A)) (5.4) 



dg \g l-gB2{X)J"'' dX 1 - g B^iX) 
and B^'{X) and B2'{X) are the derivatives of -Bi(A) and B2{X) with respect to their argu- 



ments: 
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B,'W = li^^y-'^' (2rf id + l)$^/2_i(-2A) - 24d {d - l)$^/2(-2A)) (5.5) 

B,'{X) = - i (47r)i-^/2 {2d {d + 1) $^/2_i(-2A) - 24d {d - l)%,{-2\)) 

Since these equations make no use of the fixed point values A* and g* they can be used 
to investigate both the trivial and the non-trivial fixed point. The eigenvalues and right 
eigenvectors of B, evaluated at the corresponding fixed point, then determine its critical 
exponents and scaling fields, respectively. Since, generically, B is not symmetric, its 
eigenvalues are not real and the eigenvectors are not orthogonal in general. We define the 
stabihty coefficients 6'\ 1=1,2, as the negative eigenvalues of B satisfying the equation 
BV^^ = —6^V^, where are the right eigenvectors of B. 



2. The Trivial Fixed Point 



Substituting A* = 0,(7* = into ( ^.2|) the stability matrix simplifies to 



B 



GFP 



-2 (47r)i-'^/2d(rf-3)<|.i/,(0) 
d-2 



(5.6) 



Diagonalizing ( ^.6|) then leads to two real stability coefficients with their corresponding 
right eigenvectors: 



di = 2 with 



02 = 2 -d with 



V 



(47r)i-'^/nd-3)$i/,(0) 
1 



(5.7) 



The results ( ^.7|) can be used to write down the linearized renormalization group fiow 
of the coupling constants Xk,gk in the vicinity of the Gaussian fixed point: 
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Type sign of ai asymptotic behavior 

la oi < ai-contribution dominates in 

Afc — oo in the hmit A; ^ 



Ila ai = a2-contribution dominates in (|5^ 



Afc is proportional to A;*^ ^: ^ in the hmit A; — > 0. 



Ilia ai > «! -contribution dominates in (5 

Afc — > +00 in the hmit A; ^ 0. 

TABLE I. Renormahzation group flow of Afc in the vicinity of the trivial flxed point, de- 
pending on the sign of ai. 

Xk = a^ — + a2 {Any-''/' {d - 3) <^l/M + ■■■ (5.8) 

9k = «2 + ■■■ (5.9) 

Here ai and 02 are constants of integration allowing to adjust the solution to given initial 
conditions. 

The equations ( |5.(j| ) and (§]^) show that, in c? = 4, the V^^-direction of the Gaussian 



fixed point is attractive for A; — * 0, leading to a vanishing gk, while the \^^-direction is 
repulsive. The behavior of Xk crucially depends on the sign of ai. For ai < 0, ai > 
and ai = 0, the trajectories start to the left, to the right, or on the ^^-axis. They will 
be referred to as trajectories of Type la, Ila and Illa, respectively. The corresponding 
renormahzation group behavior of A^ is summarized in TABLE ^ 

For the dimensionful coupling constants the equations (|5.8|) and ([5.9|) read 



Gk — Go + ■ ■ ■ 

A, = Ao + (47r)i-'^/2 _ 3) ^i^^^q^ Gq fc'^ + • • • (5.10) 

Here we chose M = mpi by setting 0:2 = 1, and we identified Aq = «i '^pi- From ( p.lOp it 
is easy to see that both Gk and A^ run towards constant but non-zero values in the limit 
A; — > 0, unless we set ai = by hand. But since there is no compelling reason for ai to 
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be zero we see that Aq depends on the free parameter ai and therefore on the trajectory 
chosen. Hence the Gaussian fixed point does not determine the value of the cosmological 
constant Aq in the infrared. 



3. The Non-trivial Fixed Point 

The existence of a non-trivial fixed point with g* ^ implies d — 2 + rii^{X*,g*) = 
in order for f3g to vanish. Using equation ( p.9|) this relation can be solved for g* as a 
function of A*: 

= {d - 2) B^iX*) - B,{X*) ^^-^^^ 

This equation can be used to eliminate the (yf-dependence of (3x at the non-trivial fixed 
point, i.e. we have to solve f3 ^^{X* , g* {X*)) = for A*. Due to its complicated structure 
this equation can only be solved numerically. For the /3-functions with sharp cutoff ( |3.15| ) 
the numerical evaluation yields a non-Gaussian fixed point at 

A* = 0.330, c/* = 0.403 (5.12) 



In order to discuss its properties, we take the general stability matrix ( |5.3| ) and sub- 
stitute the condition for the non-trivial fixed point, riN{X*,g*) = 2 — d. Using the sharp 
cutoff we find the following matrix entries: 



dX 



" ^ , d{d+l) \ g* 

NGFP" ^[ 2 r{d/2 + l)) l-g*B,{X* 



(4 



TT 



l-d/2 



2d{d+l) 1 12d{d-l) 1 



T{d/2-l) 1-2X* T{d/2) (1 - 2A 



*^2 



l-d/2* 2d{d+l) 



T{d/2) 1-2X* 



dg 



sc 



NGFP 



. M d(d+l) \/l B2(X*Y'' 
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2 V r(c//2) ^ ' ^ ' r(rf/2 + i) 

a/5, 



_ (47r) 1-^^/2 g*" /2d{d+l) 1 12d(rf-l) 1 \ 

NGFP ~ 3 1-g* 52(A)- l^r(d/2-l) 1-2A* r(d/2) (l-2A*)2j 

-P-^'^^ /'^taU (5.13) 

NGFP V I - g* B2[\*Y'' ) 



i?i(A)- and i?2(A)'"^ are given by equation ( |3.17| ). Substituting the numerical values of the 



non-trivial fixed point (|5.12|) and diagonalizing the stability matrix leads to the following 
pair of complex stability coefficients: 

e"- = e' + ie"= 1.941 + ^3.147 

e^ = e'- i9"= 1.941 - i 3.147 (5.14) 

In analogy to the equations (|5.8| ) and ( pl9|) the solution to the linearized fiow equation 
in the vicinity of the non-trivial fixed point then reads 

/ 



tti sin(-0" t) e"^'* Re(V^) + 02 cos(-^" t) e"^'* Im(y^) (5.15) 



9{t)-g* 

\ m - A* 

Here one sees that due to the positive real-part of 9^ the non-trivial fixed point is UV 
attractive in both directions of the (7-A— plane, i.e. attractive for t = \n{k/k) 00. 
Furthermore, the non-zero imaginary part of the stability coefficient causes the trajectories 
to spiral into the non-Gaussian fixed point when t 00. 

We emphasize that the nonzero imaginary part of the stability coefficients is not an 
artifact of our singular cutoff but appears for the smooth exponential cutoff as well, as is 
shown in TABLE |l[ 

Let us now turn to the investigation of the cutoff scheme dependence of the results 
obtained above. A change of the cutoff function R^^'^ generally leads to a change in the 
RG fiow. By definition, universal quantities, i.e. quantities that should in principle be 
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FIG. 9. The location of the non-trivial fixed point depending on the shape parameter s. While A* 
and g* show a fairly strong cutoff scheme dependence, their product X*g* is constant with a remarkable 
precision, approaching the same "plateau" values of X*g* for large values of s. In the third diagram A* is 
plotted against g* . Both of the two lines start at s = 0.8 and end at s = 15. From the different lengths of 
the lines for the sharp and the exponential cutoff one sees that the sharp cutoff leads to a much smaller 
s-dependence than its exponential counterpart. 

observable, are cutoff sclieme independent in an exact treatment. However, an artificial 
scheme dependence can arise due to the approximations one has to make in all practical 
calculations. Analyzing the cutoff dependence of universal quantities therefore provides a 
useful tool for judging the quality of the truncation. One expects that a truncation which 
yields a good approximation of the exact RG flow leads to universal quantities which are 
fairly independent of the cutoff scheme used. 

Natural candidates for testing the quality of the Einstein-Hilbert truncation are the 



stability coefficients 6' and 6" [p2|. Furthermore, one can argue that the product X*g* 
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should be universal, too |5|,|23||. The argument is as follows: The functions k ^ g^^X^ and 
their UV limits g* and A* depend on the cutoff operator and are not directly observable 
therefore. While k and, as a consequence, Gk and at a prescribed value of k cannot be 
measured separately, we may invert the function k ^ Gk and insert the result k = k{G) 
into Afc. This leads to a relationship between the Newton constant and the cosmo logical 
constant which, at least in principle, could be tested experimentally: A = X{G). In general 
this relation depends on the renormalization group trajectory chosen. But in the fixed 
point regime all trajectories approach A^ = A*/c^ and Gk = g*/k'^ which gives rise to 
\{G) = g*\*/G for A ^ mpj and G <^ rrip^. Assuming that A and G have the status 
of observable quantities, this relation shows that g*X* should be observable, and hence 
i?fc-independent. 

Therefore we also include the quantity A*^'* when we plot the numerical values for g* 
and A* obtained from the exponential and the sharp cutoff. The results are shown in FIG. 

i 

Here one finds that the non-trivial fixed point exists in the entire region under inves- 
tigation, 0.8 < s < 15. While the numerical values for A*(s) and g*{s) vastly differ for 
the various values of s and between the two families of cutoff functions one finds that the 
product X*{s)g*{s) is almost s-independent and that its values for both types of cutoff 
functions are very similar. For the "plateau" values of X*g* we find with the different 
functions R^^^ employed in this paper: 

A*^* = 0.14 for TYPE A, sc, a = 1, X*g* = 0.14 for TYPE A, Exp, a = l (5.16) 

The analogous results found in ^ using the TYPE B cutoff scheme are: 

A*^* = 0.12 for TYPE B, Exp, a = 1, A*^* = 0.14 for TYPE B, Exp, a = (5.17) 

We see that these values are very close to our present findings. In particular, the differences 
between the TYPE A and TYPE B cutoff is comparable to the difference between the 
gauge fixing parameters a = and a = 1 This supports the assumption that 
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sharp cutoff 








exponential cutoff 




s 


A* 


q* 


A*g* 


6' 


6" 


A* 


9* 


X*g* 


9' 


9" 


0.8 


0.340 


0.378 


0.129 


2.141 


3.438 


0.390 


0.233 


0.091 


1.376 


4.710 


1 


0.330 


0.403 


0.133 


1.941 


3.147 


0.359 


0.272 


0.098 


1.422 


4.307 


5 


0.281 


0.507 


0.143 


1.348 


2.743 


0.154 


0.834 


0.128 


1.499 


3.224 


10 


0.273 


0.521 


0.142 


1.294 


2.654 


0.098 


1.378 


0.144 


1.518 


3.111 


30 


0.268 


0.529 


0.142 


1.270 


2.592 


0.044 


3.149 


0.140 


1.562 


3.064 



TABLE II. Comparison between tlic numerical results for the non-trivial fixed point using the 
sharp and exponential cutoff with different values for the shape parameter s. 



the product X*g* has a universal meaning. The rather small numerical deviations of the 
plateau height can be attributed to using only a very simple truncation. 

In order to round up our survey of the properties of the non-trivial fixed point in 
d = 4, we now treat (pi and as two independent, positive numbers obeying no further 
constraints. The resulting numerical values for X*,g*,6' and 6" are shown in FIG. 0. 

This analysis clearly shows the existence of the non-trivial fixed point for the complete 
region under investigation, < (pi,(p2 < 10. Even though the magnitude of the stability 
coefficients varies by about a factor of three, as an effect of our truncation, the qualitative 
properties of the fixed point are the same for all values of ipi and (p2- The fixed point always 
possesses a pair of complex-conjugate stability coefficients whose real-part is positive, 
Re(6'^) = 9' > 0, i.e. the fixed point is always UV-attractive. This provides us with 
further evidence for the existence of a non-trivial fixed point in the full theory 

B. The Phase Portrait 

Motivated by the very good agreement between the exponential and the sharp cutoff 
found when solving the M-scaled flow equation in Section III and analyzing the non-trivial 
fixed point, we now investigate the full A-f?— parameter space, using the /3-functions (|3.15|) 
with the technically much more convenient sharp cutoff. 
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FIG. 10. Numerical results for the location and stability coefficients of the non-trivial fixed point, 
treating ipi and ip2 as independent. The solid line corresponds to the results obtained by considering ipi 
and (p2 as a function of the shape parameter s. For all values of ipi and ip2 the non-Gaussian fixed point 
is seen to have the same qualitative features. 

But before presenting the numerical solutions, we summarize the properties of the 
parameter space that can be directly read off from (|3.15|) . As has been noted in the 
context of ( [4. 81 ), the flow equation is singular for \ > y/2. In terms of the /c-scaled 
coupling constant = \{y)/y this singularity occurs for A > 1/2, resulting in a boundary 
of the A-fyf— parameter space given by the line A = 1/2. 

The vanishing of /3g{X, g) for g = and arbitrary A leads to a separation of the phase 
space into two decoupled regions with positive and negative coupling g, respectively. 
Trajectories starting in one of these regions will never cross the separation line g = 0. 

The characteristics of the trivial and the non-trivial fixed point have already been 
discussed in the previous subsection. One expects that the non-Gaussian fixed point 
dominates the RG flow for large t = ln{k/k), leading to a spiraling in of the trajectories 
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FIG. 11. The points in the parameter space where jy^r diverges for the sharp and exponential cut- 
off. The bars indicate the location of the non-trivial fixed point for the exponential cutoff with shape 
parameters s — 1 and s = 30. 

on the point {X*,g*). The stabihty axis V"^ of the Gaussian fixed point will provide 
a separation between the trajectories which run towards = — oo and = +1/2, 
respectively. 

In terms of the coordinates X,g it is also easy to visualize the singularity of the 
anomalous dimension rj^ encountered in Section III. The function riN{X,g) diverges at 



those points {g, A) at which the denominator of equation (|2.9| ) vanishes. The function 



g{X) defined by 1 — g{X)B2{X) = is shown in FIG. |TT]. For the sharp cutoff one finds 
that this singularity occurs on a s-independent horizontal line at 5^ = — Gvr/S. For the 
exponential cutoff the divergence of 7]^^ leads to a (s-dependent) line in the first quadrant 
of the X-g-plane. 

We shall now derive the full phase portrait of the Einstein-Hilbert truncation by 
numerically solving the autonomous differential equation for the "fc-scaled" quantities gk 
and Afc. In the first step we focus on the part of the A-t^-plane where the RG fiow is 
governed by the two fixed points found in the previous subsection. 



The resulting fiow diagram is shown in FIG. |12. 

This figure clearly shows the separation between the trajectories with positive and 
negative coupling g. 

The most interesting feature of this diagram is the interplay between the non-Gaussian 
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FIG. 12. Part of the parameter space with its RG flow. The arrows along the trajectories 
point in the direction of the renormalization group flow, i.e. towards decreasing values of k. The 
flow pattern is dominated by a non-Gaussian fixed point in the first quadrant and a trivial one 
at the origin. 



and the Gaussian fixed point in the positive coupling region. In the limit t —>■ oo the 
non-Gaussian fixed point completely dominates the flow of all trajectories in this region, 
leading to their spiraling in on the point (A*, g*). This result exactly matches the behavior 
expected from the complex stability coefficients found in the previous subsection. The 
behavior of the trajectories in the vicinity of the trivial fixed point confirms TABLE 
|: Trajectories starting to the left of the separation line (ai < 0) are running towards 
A —>■ —oo while those to the right (ai > 0) terminate at the border singularity A = 1/2 
at finite values of k. 

The trajectory separating these two regions is of special interest and will be called 
"the" separatrix. It connects the non-trivial fixed point in the UV to the trivial fixed 
point for /c — s> 0. As a result, this trajectory leads to a vanishing renormalized cosmological 
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FIG. 13. Full phase structure of the Einstein-Hilbert truncated theory. On the bold hori- 
zontal line ry^ diverges. 

constant Aq- 

The trajectories running to the left of the separatrix also possess a well defined limit 
k 0. They lead to negative values of the renormalized cosmological constant Aq- 

In the next step we now extend our numerical survey to the full parameter space. The 
resulting renormalization group flow is shown in FIG. 0. Here it becomes obvious that the 
negative-(7 region also contains a trajectory which separates the regions with trajectories 
running towards A —>■ —oo and A 1/2. Furthermore one sees that the singularity of 

at (7 = — 67r/5 results in a separation between trajectories showing a screening and 
ant i- screening behavior of the Newton constant in the IR. This resembles the behavior 
found for the exponential cutoff in Section III. 



In a first step of classifying the trajectories found in FIG. 13, we project the renormal 



ization group trajectories of the full system onto the g-axis. The results are displayed in 
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FIG. 14. Projecting the full renormalization group flow onto the g-axis leads to a separation between 
the trajectories in the strong, weak and negative coupling region. 



FIG. |T^. According to their limit for k ^ or t = \n{k/k) —>■ — oo three different classes 
of trajectories can be distinguished: 

(i) Trajectories with lim g{t) oo. They form the "strong coupling region". 

t — ^— oo 

(ii) Trajectories with lim g(t) = 0. They form the "weak coupling region". 

t — oo 

(iii) Trajectories with lim g{t) < 0. They form the "negative coupling region". 

t — oo 

Note the oscillating behavior of g{t) before the trajectories adopt their asymptotic value. 
This is caused by the nonzero imaginary part of the stability coefficients found for the 
non-trivial fixed pointQ. 



In a second step we classify the trajectories shown in FIG. |T3| according to their starting 
and end points. FIG. |1^ shows the resulting phase space regions, which are distinguished 
by a different kind of dashing of their trajectories. 



*A plot of g{t) similar to FIG. 14 has been given in ref. see fig. 2 there. In this reference 
the running of A has been disregarded, however. As a consequence, no oscillations were found. 
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-0.4 -0.2 0.2 0.4 

FIG. 15. Classification of the full Einstein-Hilbert truncated theory space. The trajectories Ila and 
lib end at the trivial fixed point and separate the regions la and Ilia, lb and Illb, respectively. They are 
not labeled explicitly. 

The characteristics of each region are summarized in TABLE ^TT| which contains the 
classification of all trajectories occurring in the Einstein-Hilbert truncation. TABLE 
pni is organized as follows: The first column labels the type of the trajectory as it is 
marked in the phase space diagram FIG. [1^. (Only the single trajectories of Type 11a 
and lib separating the regions la and 111a, and lb and Illb, respectively, are not marked 
explicitly in this diagram.) The columns labeled "UV-" and "IR-behavior" indicate the 
characteristic features of the trajectories, "UV" referring to the endpoint of the trajectories 
for k oo and "IR" relating to 0. These limits do not exist for all the classes. The 
aborting of the trajectory at a finite value of k in either the UV or the IR is indicated 
by "Sing.". The values of X,g given in the table indicate where the corresponding RG 
trajectories end. 
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TABLE III. Classification of all RG trajectories occurring in the Einstein-Hilbert truncation. 



Type 



UV-behavior 



IR-behavior coupling region 



la 



Ila 
separatrix 



Ilia 



IVa 



Va 



lb 



lib 



Illb 



NGFP {X*,g*) o.k. 



NGFP {X*,g*) 

NGFP {X*,g*) 

NGFP {X*,g*) 
X = 1/2 



Sing. 



Sing. 



Sing. 



Sing. 



g>0 

X = 1/2 
g<0 

X = 1/2 
g<0 

A = 1/2 
g<0 



GFP 



Xk -oo 
9k ^0+ 



Xk^O+ 
9k ^0+ 



Sing. 



Sing. 



Sing. 



o.k. 



GFP 



Sing. 



A = 1/2 
g>0 



Xk -oo 
9k ^ oo 



Xk -oo 
9k ^ oo 



gk^O 
Xk — — oo 



9k 0" 



A = 1/2 
g<0 



weak coupling 



weak coupling 



weak coupling 



strong coupling 



strong coupling 



negative coupling 



negative coupling 



negative coupling 



IVb 



Vb 



Sing. 



Sing. 



g = — 67r/5 



A 

g = — 67r/5 



o.k. 



Xk — — oo 
9k 0" 



Sing. 



A = 1/2 
g < — 67r/5 



negative coupling 



negative coupling 
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In the column "UV-behavior" the label "NGFP" means that the trajectory runs into 
the non-Gaussian fixed point. The trajectories labeled with "Sing." either end at the 
boundary line A = 1/2 where we distinguish between the two regions g > and g < or 
in the singular line caused by the divergence of rf^ dX g = — Gvr/S and arbitrary values of 
A. Note that the limit k oo exists only for the trajectories running into the "NGFP". 

In the IR, trajectories which possess an IR-limit k are indicated by the note 
"o.k.". They lead to negative values of the renormalized cosmological constant Aq with 
positive {gk 0+) or negative {gk 0^) Newton constant, respectively. The label 
"GFP" indicates that the corresponding trajectory ends at the Gaussian fixed point, 
yielding a vanishing cosmological constant Aq. Singular behavior, i.e. the termination of 
the trajectory at a finite value of k, appears when the trajectory reaches the boundary 
A = 1/2 in the positive {g > 0) or negative {0 > g > — Gvr/S, g < — Gvr/S) coupling region. 
The possibility of gk diverging at a finite value of k is indicated hj gk ^ oo, Xk ^ —oo. 

The column "coupling region" finally specifies to which coupling region, as defined in 



FIG. the trajectory belongs. 

The Newtonian regime consists of the trajectories in the weak coupling region. They 
are the only ones yielding a finite, positive Newton constant Gk in the IR. In this respect 
it is important to note that only the trajectories of the Types la and Ila can be continued 
down to = 0, leading to a finite value of Go and a negative or vanishing cosmological 
constant, respectively. Since for these two classes the limit A; ^ oo is also well-defined, 
these solutions could possibly be used in order to define a fundamental quantum theory of 
gravity with a vanishing or a negative renormalized cosmological constant, respectively. 
The trajectories of Type Ilia, which run towards positive values of A, terminate at the 
singularity A = 1/2 and do not give rise to well-defined renormalized parameters Aq and 
Go at k = 0. 
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C. Crossover Behavior 



The phase diagram FIG. |T2| shows that the trajectories la, Ila and Ilia which are 
relevant for the Newtonian limit of quantum gravity cross over from the non-trivial fixed 
point in the UV to the basin of attraction of the trivial fixed point in the IR. In this 
subsection we first concentrate on the classes la and Ila. They lead to finite A; ^ limits 
of the coupling constants. Since these trajectories extend down to /c = they allow us to 

— 1/2 

identify the scale M with the Planck-scale mpi = Gq . 

Due to the trivial and the non-trivial fixed point governing the RG flow of these 
trajectories in the limits k —>■ and k —>■ oo, respectively, the trajectories of these classes 
obey simple scaling laws in these regimes. 

In the UV one finds the following scale dependence of the coupling constants, valid 
for both types of trajectories: Afc const, gk const -v^ Gk, G oc A^, A oc k"^. 

In the IR the trivial fixed point leads to the scaling laws given by the equations (|5.8|) 
and (p.9|). For Type la trajectories they are ex. k"^ ^ Gk, G oc oc /c"^ -v^ A/,., A oc 

/c°, and for the Type Ila they read gk oc k'^ Gk, G oc k^, A^ oc A;^ -v^ \k, A oc k^. 

Next we determine the scale k, at which the transition between these regimes takes 
place. Thanks to the identification M = mpi the related dimensionful quantity k = kmpi 
has a clear physical interpretation. We solve equation ( p.l5| ) with the "final condition" 
G{k = 0) = 1, together with A(A; = 0) = —0.01 and A(A; = 0) = for a typical trajectory 



of Type la and Ila, respectively. This leads to the results displayed in FIGs. ^ and |T7. 

The main conclusion is that the change in the scaling laws arising from the two fixed 
point regions is located rather close to the Planck-scale, k = mpi. For k ^ 0.1 and k ^ 10 
the respective scaling laws of the trivial and non-trivial fixed point already dominate the 
running of the coupling constants. 

This picture is confirmed by the diagrams displaying the anomalous dimension k 



ri]\f{Xk, gk) along the trajectories la, Ila and Ilia shown in FIG. |T^. In the region governed 
by the non-trivial fixed point one has rji^j ^ —2 while in the IR-region T]j\f vanishes. 
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PI 



FIG. 16. Scaling laws for Typ la trajectories. For k < 0.1 and fc > 10 the flow is governed by 
the scaling laws of the Gaussian and the non-Gaussian fixed point, respectively. The transition between 
these scaling laws happens at A; w mpi. The cusp appearing in the double logarithmic diagrams for the 
modulus of the cosmological constant is caused by X{k) and Afe becoming negative at finite k. 
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log Jg(k)) log. JX(k)) 




FIG. 17. Scaling laws for the trajectory Ila, the separatrix. For k < 0.1 and > 10 the flow is 
governed by the scaling laws of the Gaussian and the non-Gaussian fixed point, respectively. The change 
between these power laws happens at A; w mpi. 
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FIG. 18. Anomalous dimension yyjv along the Type la, Ila and Ilia trajectories shown in the first 
diagram. 

(Recall that 7]*^ = —2 and r/]^ = at the non-Gaussian and the Gaussian fixed point in 
d = 4, respectively.) For the Type Ilia the termination of the trajectory is accompanied 
by a steep decrease of 1]^, caused by the divergence of Bi{XY^ at constant B2{Xy'^ for 
A approaching the boundary line A = 1/2. (This can easily be checked from equation 
( ^.17] ).) The vast decrease of 1]^ for the trajectories Ilia suggests that the Einstein- 
Hilbert truncation may not be sufficient to describe the RG fiow close to the boundary 
line A = 1/2. 

Motivated by the important role played by the product \*g* we also plot k i-^ \{k)g(k) 
for those trajectories which were considered in FIG. |l^. The results are shown in FIG. 
EOl. 
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FIG. 19. Product X{k)g{k) for the representative trajectories la, Ila and Ilia shown in FIG. |l^. The 
parameter k is given by A; = fc/mpi for Type la and Ila, and k = k/M for Type Ilia. 

Here we observe that for small values of fc, A; ^ 1 say, the value of \{k)g(k) along the 
separatrix provides a clear separation of the trajectories la and llla. In the region above 
k ^ 1 the trajectories cross. In this region the Af^-value along the separatrix does not 
provide a separator between the trajectories la and llla. This behavior is caused by the 
non-zero imaginary part of the stability coefficients of the non-Gaussian fixed point. It 
leads to the observed oscillations of the trajectories in the intermediate region. In this 
respect we point out that for a non-Gaussian fixed point which by chance has lm(^^^) = 0, 
the quantity \{k)g{k) along the separatrix would provide a good separator between the 
trajectories running in the regions la and llla for all values of k. This remark will become 
important when we compare our results to those found in lattice calculations in Section 
VI. 

D. The Non- Trivial Fixed Point in other Dimensions 

The /3-functions ( |2.(j| ) and the entries of the stability matrix B in ( p. 3D have been 
derived for arbitrary d and allow for an investigation of the non-trivial fixed point in 
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FIG. 20. Comparison of the numerical values of X*,g*,0' and 6" for different cutoff functions in 
dependence of the dimension d. The upper line shows that for 2 + e < d < 4 the cutoff scheme dependence 
of the results is rather small. The lower diagram shows that increasing d beyond about 5 leads to a 
significant difference in the results for 9' , 6" obtained with the different cutoff schemes. 

any dimension d. For 2 + e < c? < 4, the existence of the non-trivial fixed point for the 
/3-function with the exponential cutoff (|3.2|) has already been demonstrated in 

FIG. ^ shows the numerical values of X*,g*, 9' and 9" obtained by using both the 
exponential and the sharp cutoff. The s-dependence of the sharp cutoff is again defined 
by the relation ( p.l3|) which trivially extends to any d. 

The results in FIG. ^ show that the scheme dependence of both the (nonuniver- 
sal) location and the (universal) critical indices at the fixed point increase steadily with 
increasing d. But up to dimensionalities between about 4.5 or 5, say, one still finds a rel- 
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FIG. 21. Numerical values A* and (g*)^^ for the non-trivial fixed point in d > 4. The branching 
occurring for the sharp cutoff indicates the existence of a second fixed point in the negative coupling 
region. 

atively good agreement between the results obtained with the different cutoff functions. 
Beyond this point the universal quantities 9' and 6" become seriously scheme dependent. 

It is reassuring to see that there seem to be no qualitative differences between d = A 
and d = 2 + e. We interpret this as another confirmation of the reliability of the Einstein- 
Hilbert truncation in 4 dimensions, at least at a qualitative level. 

Looking at the cutoff dependence of A*, (7* in the region (i ^ 4, one finds that the 
trend of the growing scheme dependence of the fixed point properties continues. FIG. ^ 
shows that the value of g* vastly increases for increasing d, and that A* approaches the 
singularity at A = 1/2. Once A* = 1/2 is reached, the non-trivial fixed point disappears 
at a certain critical dimension dc which strongly depends on the cutoff function. In FIG. 
2TI one finds: 



(ic ~ 6 for the exponential cutoff with s = 1. 
(ic ~ 17 for the sharp cutoff with s = 1. 



dc > 26 for the sharp cutoff with s = 30. 



(5.18) 



The first diagram of FIG. ^ shows quite impressively that d = 4 seems still to lie 
on the safe side of a rather pronounced "phase-transition" at ^ 5. At this point, g* 
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FIG. 22. f3)^{\, g* {X)) as a function of A for the exponential and the sharp cutoff with shape pa- 
rameter s = 1 for selected dimensions d. The zero of this function is A*. The left diagram demonstrates 
the escape of the non-trivial fixed point into the boundary A = 1/2 for the exponential cutoff. The right 
diagram shows the appearance of a second zero and the vanishing of the familiar non-trivial fixed point 
for the sharp cutoff. 



suddenly jumps from g* < 1 to g* ^ 1, so that we must expect the truncation to become 
problematic at d ^ 5. 

In order to better understand the disappearance of the non-Gaussian fixed point, we 
plot in FIG. the function A ^ Px{^, g*{^)), whose zero is A*, for both the exponential 
and sharp cutoff and for various dimensions d. The function g*{X) thereby is given by 
equation ( ^.111 ). One finds that, for the exponential cutoff, f3^{X, g*{\))^^^ is positive 
definite for d > d^. which implies that there can be no non-trivial fixed point. For the 
sharp cutoff, /3a (-^j fi'*('^)) develops a second zero in the negative coupling region before the 
old non-trivial fixed point in the positive region escapes through the A = 1/2 boundary. 
The appearance of the second zero is indicated by the branching shown in the \*-d- 
diagram of FIG. |2^. For d > dc only the new fixed point in the negative coupling region 
remains. 

The properties of these fixed points for selected dimensions are shown in TABLE |V| 
and TABLE for the sharp and the exponential cutoff, respectively. 

The data shown in TABLE and suggests that the Einstein-Hilbert truncation 
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TABLE IV. Numerical values for the non-trivial fixed point and the fixed point appearing in the 
negative coupling region for the /3-functions with sharp cutoff with s = 1 in selected dimensions d. Here 

'-' denotes that no fixed point of the corresponding type exists for this dimension. A pair of real- valued 
stability coefficients, appearing for c? > 6, is indicated by giving the two corresponding O-YaXiica in the 
^'-column. For d > 6 the results found here are certainly not reflecting properties of the full theory 
and should be seen as a demonstration of the limitations of the Einstein-Hilbert truncation in higher 
dimensional spacetimes. 







non- 


trivial fixed point for exponential cutoff with s 


= 1 
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A* 


9* 


\*g* 


e' 


e" 
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0.140 


0.133 
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1.109 
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0.359 


0.272 


0.098 
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4.307 
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0.445 


0.636 


0.283 


3.419 


8.503 
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0.497 


1.673 


0.831 


10.02 


16.82 



TABLE V. Numerical values for the non-trivial fixed point with exponential cutoff in selected 



dimensions d. For d > 6 no non-trivial fixed point exists. 



46 



should produce reliable results up to c? = 4 or maybe d = 5. Beyond that the results shown 
in TABLE ^ are untrustworthy and should be understood as an illustration of the effects 
arising from an improper truncation. For the non-trivial fixed point this unreliability is 
indicated by the fixed point lying very close to the boundary line A = 1/2 and at large 
values g* ^ 1. In this region of the A-t?— plane, the termination of the trajectories and the 
steep decrease in r^^r indicate that the Einstein-Hilbert truncation might be insufficient in 
describing the RG flow. The appearance of the second zero of the /3-function with sharp 
cutoff and the absence of this new zero for the smooth exponential cutoff point towards a 
strong scheme dependence even at the qualitative level. This is a typical symptom of an 
insufficient truncation. 

Summarizing the results of this subsection we find that the Einstein-Hilbert truncation 
is most likely insufficient to describe the RG flow for d ^ 5. This limitation is due to 
operators like / d'^x^/gR'^ and higher powers of the curvature scalar which are not included 
in the truncation. Based on the canonical dimensions of their coupling constants one 
expects an increasing importance of these terms in the UV as d is increased. Therefore 
it is likely that in order to properly describe quantum gravity in higher dimensional 
spacetimes a more refined truncation will be needed. 

VI. DISCUSSION AND CONCLUSION 

In this paper we studied the exact renormalization group equation for the effective 
average action of pure quantum gravity in its original formulation [Q]. The RG flow 
of Newton's constant and the cosmological constant was investigated nonperturbatively 
within the Einstein-Hilbert truncation of the theory space. The new contributions of 
the present paper are the introduction of a sharp cutoff and a comprehensive numerical 
analysis of the RG trajectories in 4 spacetime dimensions. 

As in ref. 0, we used a cutoff action A^S whose general structure is of "TYPE A" 
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in the terminology of ref. ||^[]]. It contains a shape function R^^^ for which we considered 
both a smooth exponential ansatz and, as a singular limit, a sharp step function of 
infinite height. While in other applications sharp cutoffs often lead to inaccurate or 
even undefined or divergent results we found that the flow equations for g and A remain 
perfectly well defined in the sharp cutoff limit. Moreover, all RG trajectories which we 
computed with both the smooth and the sharp cutoff turned out virtually identical. Thus 
having confidence in the technically much more convenient sharp cutoff we employed it for 
a complete classification and computation of the RG trajectories on the (7- A— parameter 
space. Our main results are summarized in FIGs. |1^ and |I^ and in TABLE |T|. 

The most prominent feature of the RG flow resulting from the Einstein-Hilbert trun- 
cation is a non-Gaussian fixed point which acts as an UV-attractor for the trajectories la, 
Ila, Ilia and IVa on the g > half plane. It is an extremely important question whether 
this fixed point is a truncation artifact or whether it is also present in the full theory. In 
the latter case 4-dimensional Quantum Einstein Gravity is likely to be "asymptotically 
safe" [^], i.e. nonperturbatively renormalizable by taking the infinite cutoff limit at this 
fixed point. Einstein gravity would then have the status of a fundamental rather than 
merely an effective theory, and it could be valid at arbitrarily small distance scales. 

In refs. [^,0 this question was investigated in detail. The reliability of the Einstein- 
Hilbert truncation was tested both by analyzing the scheme dependence within this trun- 
cation 1^ and by adding a higher-derivative invariant to it [0. The picture suggested 
by these investigations is that, in 4 dimensions, the RG flow in the vicinity of the non- 
Gaussian fixed point is very well described by the Einstein-Hilbert truncation. Hence we 
have very good reasons to beheve that the fixed point actually should exist in the full 



^This cutoff should not be confused with the new one introduced in ref. whose general struc- 
ture is of the "TYPE B" which appears naturally when the transverse-traceless decomposition 
of the metric is used. 
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theory. Only when one lowers (!) the cutoff and leaves the region of asymptotic scaling 
more complicated operators such as R^, R^, ■ ■ ■ are generated. 

The results of the present paper provide further evidence supporting the hypothesis 
of the non-Gaussian fixed point. In order to judge the reliability of the truncation, we 
checked the cutoff scheme dependence of various universal quantities {6',6",g*X*) which 
are expected to be scheme independent in an exact treatment. Quite remarkably, the 
results found here for a sharp cutoff of TYPE A are rather similar to those found in |^ 
with a smooth cutoff of TYPE B. Typically the scheme dependence within the family 
of sharp cutoffs (parameterized by ipi and ip2) is of about the same magnitude as the 
differences between the sharp and the smooth cutoff, and between the TYPE A and 
TYPE B structure. In particular the product X*g* was found to be scheme independent 
with a very surprising precision. Here the sharp cutoff leads to even slightly more stable 
results than the smooth one, as is shown in the last diagram of FIG. |^. 

We also investigated how the scheme dependence of the fixed point data varies with 
the dimensionality d. With increasing d the quality of the Einstein-Hilbert truncation 
deteriorates and it certainly becomes insufficient at about d = 6 where the existence or 
nonexistence of the non-trivial fixed point depends on the cutoff chosen. It seems that 
in (i = 4, at least qualitatively, the conditions are still very similar to d = 2 + e with 
< e ^ 1. All admissible cutoffs lead to the presence of the non-trivial fixed point. They 
agree on a quantitative {X*g*) or at least semi- quantitative level {6', 6"). 

While in [^,^ only the linearized flow near the fixed points is discussed, the numerical 
investigations of the present paper allow us to follow the trajectories emanating from the 
non-Gaussian fixed point all their way down from the UV to the IR. In particular we were 
able to determine the scale fcasym where the asymptotic scaling behavior governed by the 
UV fixed point comes to an end. For the trajectories which can be continued to A; = we 
found that /casym ~ ""^pi = G{k = 0)~^/^. Below this scale the RG flow is governed by a 
trivial fixed point at the origin of the A-f?— plane. The (Planck) scale /casym is analogous to 
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the mass scale Aqcd in QCD in the sense that it marks the lower border of the asymptotic 
scaling region. 

Probably the analogy between QCD and gravity goes even further. For k ^ Aqcd, 
thanks to asymptotic freedom, the effective average action of Yang-Mills theory has a 
very simple local structure, and simple truncations such as Tk = 1/4 Zk J d'^x{FI^^)^ are 
sufficient |^^. Only when k is close to, or much smaller than Aqcd its structure becomes 
very complicated and nonlocal. Also gravity seems to be asymptotically free 0, and if 
the results of and the present paper point in the right direction, gravity, too, can 
be described by a simple local truncation above /casym ~ fnpi- Only when we lower k 
down to the Planck scale, F^, will contain higher order local invariants {K^, ■ ■ ■) as well 
as nonlocal invariants [^,^]. Presumably the use of the Einstein-Hilbert truncation is 
much more problematic at ~ mp\ than at A; ^ mpi. In fact, the aborting trajectories 
of Type Ilia and the unbounded increase of |?7Ar| which precedes their termination most 
probably hint at an insufficiency of the Einstein-Hilbert truncation in the infrared. It 
must be emphasized that there are no analogous consistency problems above fcasym- 

It is likely that in the IR a proper treatment of the Type Ilia trajectories (or rather 
their counterparts in the exact theory) requires much more sophisticated truncations. It 
is plausible to speculate that those truncations will encode the strong IR-quantum effects 
which have been discussed in refs. f^. They might be the key to a dynamical resolution 
of the cosmological constant problem. 

The above picture of quantum gravity being "simple" above mpi and "complicated" 
below mpi contradicts the general prejudice that far below mpi gravity is well described 
by the "simple" Einstein-Hilbert action and becomes very "complicated" for all k ^ mpi. 
We think that this prejudice might turn out wrong for (at least) two complementary 
reasons: (i) Ordinary perturbation theory predicts that F^ becomes "complicated" (higher 
operators are generated) when k approaches the Planck scale from below. While this is 
probably correct, it is an unjustified extrapolation that continues to be "complicated" 
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for k ^ mpi- In this regime the "simphcity" due to asymptotic freedom sets in, but this 
cannot be discovered by ordinary perturbation theory, (ii) The classic experimental tests 
of general relativity all refer to the theory without a cosmological constant. Prom this 
side nothing is known about the IR-properties of the (quantum) theory with A which 
could be rather "complicated" without contradicting any experimental fact. 

The weak coupling region of the A-^^— plane contains the trajectories of Type la, leading 
to Ao < 0, and of Type Ilia, terminating at the boundary A = 1/2. The two classes of 
trajectories are separated by a single trajectory of Type Ila, the separatrix. It runs from 
the non-trivial to the trivial fixed point and leads to a vanishing renormalized cosmological 
constant Aq = at its cndpoint. 

We saw that the la and Ila trajectories can be continued down to A; = without any 
problem and that |77Ar| remains bounded along these trajectories. This might indicate 
that in this sector, characterized by a renormalized cosmological constant Aq < 0, the 
Einstein-Hilbert truncation is reliable even in the IR, at least at a qualitative level. With 
the present calculational techniques it would be very hard to check this conjecture on a 
theoretical basis. However, as far as the case Aq = is concerned, the phenomenological 
success of classical general relativity (without a cosmological constant) is an excellent 
confirmation of the Einstein-Hilbert action in the IR. 

The Einstein-Hilbert action thus being a good approximation to Tk both in the extreme 
UV and the extreme IR, it is plausible to assume that the scparatix which connects the 
two fixed points does not change much when we go from the truncated to the full theory. 
It is then possible to define the theory called "Quantum Einstein Gravity" (with zero 
renormalized cosmological constant) by means of this specific RG trajectory. 

It is natural to interpret the separatrix as a kind of phase transition line. It is in fact 
similar to the critical line running into the Gaussian fixed point in scalar theories. In the 
latter case this line separates the symmetric from the spontaneously broken phase, i.e. 
trajectories with positive and negative (mass)^, respectively. In the case of gravity the 
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line separates the trajectories which, in the IR, lead to a positive or negative cosmological 
constant, respectively. On the side with A > the RG trajectories terminate at a certain 
finite value of k. Whether this is an artifact of the truncation or a real physical effect 
cannot be decided on the basis of the present analysis. 

In this context it is intriguing to compare our results to those obtained by numerical 
simulations within Regge's simplicial formulation of gravity p7|P^ . These studies indicate 
that simplicial quantum gravity in four dimensions exhibits a phase transition in the 
bare coupling G between the following two phases: a strong coupling phase in which 
the geometry is smooth at large distance scales with (g^^) ~ c6fj,u, and a weak coupling 
phase in which the geometry is degenerate, (g^u) = (branched polymer-like phase). 
Interestingly enough, in the strong coupling phase one finds a small negative average 
curvature. 

It is tempting to identify the strong coupling phase with our Type la trajectories which 
ultimately go to negative Aq and negative curvature i? = 4A, and the weak coupling phase 
with the type Ilia trajectories which have no limit k ^ and cannot describe large, nearly 
fiat universes therefore. Hamber parameterizes the bare couplings by a quantity k 
(not to be confused with our cutoff) which reads, in our notation. 



SttGA 



-1/2 1/2 



[Svr^^A]- 



Strong (weak) coupling means k < kc (k > kc) for a certain critical value kc. This matches 
precisely with our findings in Section V. Depending on whether the product gX takes on 
negative or positive values in the IR, either the Type la or the Type Ilia trajectories are 



realized. This can be read off from FIG. 19 provided one ignores the oscillations and the 



crossing of the trajectories which are caused by the imaginary part of the critical exponent 
6. It remains to be seen if this identification is correct. 

A quantity which in principle might lend itself to a quantitative comparison is the 
exponent 6. In [^] the analogous exponent is assumed to be real, 6 = 6' = l/u, with 
u the conventionally defined critical exponent used in the theory of critical phenomena. 
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The fit to the data yields a value close to v = 1/3 so that we should expect 6' ~ 3. 
The typical values we found in Section V are somewhat smaller. However in a detailed 
investigation of the fixed point with a generalized truncation and a quantitatively probably 
more reliable smooth cutoff has been performed. The results suggest that an improved 
calculation indeed could stabilize close to 6' ~ 3. It is unclear, however, how the imaginary 
part of 6 should be properly taken into account in this comparison. 



Acknowledgement: We would like to thank A.Bonanno, W.Dittrich, O.Lauscher and 
C.Wetterich for helpful discussions. 

APPENDIX A: SINGULARITY OF rjjy FOR THE EXPONENTIAL CUTOFF 



When investigating the complete system ( p.l3| ) and ( |2.14| ) in Section IV we omitted 



discussing some peculiar properties of the RG flow derived with the exponential cutoff. 
These properties will be analyzed in this appendix. 

The key to understanding this behavior is the fact that there exists a line on the A- 
(jf— plane along which ri^^{X,g) diverges. This line is shown in FIG. O. Switching from 



the variables (\,g,k) to {X,G,y) leads to an analogous divergence of r]^^{\,G,y = k'^) 
on a 2-dimensional curved surface in the 3-dimensional A-G-y— space. We can visualize 
this surface by intersecting it with some plane G = const, which then leads to a line 
of singularities on this A-y— plane. For a given value of G, we parameterize this line as 



y ^ Asing(y). The result for G = 0.5 is shown in FIG. One finds that the line Asing(y) 
is located below the boundary A = y/2, approaching it as y — 0. This line separates 
the A-y— parameter space (for G fixed) into a region below (A < Asing(y)) and above 
(Asing(y) < A < y/2) it. In these regions one finds t]^^ < and i]^^ > 0, respectively. 
From a purely mathematical point of view, one can pose admissible initial conditions for 
the RG equations in both regions. 
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FIG. 23. Line of ?7^'^''-singularities for fixed G = 0.5. This line separates regions with > 

above it and with f]^^ < below it. 

Two typical trajectories starting below and above Asing(^) are shown in FIG. |2^. The 
most probably unphysical trajectories starting in the narrow region Xsmgiv) < ^iv) < y/2 
terminate at the boundary A = y/2. They would lead to a screening behavior of the 
Newton constant, i.e. G{y) increases with increasing y = k'^/M'^. 

The trajectories starting with X{y) < Asing(y), i-e. in the region where 7]^"^ < 0, lead to 
the RG trajectories discussed in Section IV. C which are physically relevant and show the 
expected antiscreening. According to their initial values, these trajectories are solutions 
of Type la, Ila or Ilia. 

Here it is interesting to note that trajectories starting close to the singularity of r/Tv? i-e. 
X{y) ^ Asing(y) can be continued down to "unnaturally low" values y. The corresponding 
trajectories cross the ones starting at lower values A(^) in the A-y— plane before they 
terminate at the boundary. An example is the trajectory with the largest A(y) in the first 
plot of FIG. 1^. Since all trajectories of this type terminate, this mechanism cannot be 
used to change the Type Ilia character of a trajectory. 

We want to point out that the region above the singular line at Aging (y) is probably 
unphysical. This assumption is supported by the fact that our investigation with the 




54 




y 



FIG. 24. Typical trajectories starting in the regions A < Aging (y) (dashed Hne) and 
Aging (y) < A < y/2 (sohd line). While the trajectory below Aging (y) shows the typical behavior dis- 
cussed in Section IV, the solution above Aging (y) yields a screening behavior of the Newton constant. 

sharp cutoff did not yield any comparable behavior in the positive coupling region. The 
trajectories calculated with the smooth and the sharp cutoff, respectively, agree quite 
well for y > yterm- Only close to their endpoint qualitative differences appear, such as 
the earlier termination at Aging (y) rather than A = y/2. This scheme dependence is a 
typical symptom showing that the truncation becomes unreliable close to the boundary 
of parameter space. 
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